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SUZAKU OBSERVATION OF THE BRIGHTEST BROAD-LINE RADIO GALAXY 4C 50.55 (IGR J 21247+5058) 
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ABSTRACT 

We report the results from a deep Suzaku observation of 4C 50.55 (IGR J 21247+5058), the brightest 
broad-line radio galaxy in the hard X-ray (> 10 keV) sky. The simultaneous broad band spectra over 
1-60 keV can be represented by a cut-off power law with two layers of absorption and a significant 
reflection component from cold matter with a solid angle of J7/27T « 0.2. A rapid flux rise by ~ 20% 
over 2 x 10 4 sec is detected in the 2-10 keV band. The spectral energy distribution suggests that 
there is little contribution to the total X-ray emission from jets. Applying a thermal Comptonization 
model, we find that corona is optically thick (r e ~ 3) and has a relatively low temperature (kT c w 30 
keV). The narrow iron-K emission line is consistent with a picture where the standard disk is truncated 
and/or its inner part is covered by optically thick Comptonizing corona smearing out relativistic broad 
line features. The inferred disk structure may be a common feature of accretion flows onto black holes 
that produce powerful jets. 

Subject headings: galaxies: active - galaxies: individual (4C 50.55) - X-rays: galaxies 



1. INTRODUCTION 

Radio galaxies are a class of Active Galactic Nuclei 
(AGNs) spouting powerful radio jets whose axis is not 
aligned along the line of sight. While the formation 
mechanism of relativistic jets is not fully understood, it 
must be closely related to mass accretion flow onto the 
central supermassive black hole. In fact, studies of Galac- 
tic black holes have revealed that the relative power of 
the jets to accretion critically depends on the state of the 
accretion disk, which is predo minantly determined by the 
mass accretion rate (e.g., see lFender et al.ll20 04b). For 
AGNs, however, detailed studies of accretion disk state in 
relation to the jet formation are still limited. Broad band 
observation of nearby, bright radio galaxies hence give us 
ideal opportunities to unveil this problem, because, un- 
like "blazars" , the innermost disk can be well observable 
in the X-ray band with much smaller contribution from 
the jet emission. 

4C 50.55 (IGR J21247+5058) is a Fanaroff- Riley 
type II broad line radio galaxy (BLRG), discov- 
ered by INTEGRAL in its first survey catalogue 
(|Bird. Barlow fc Bassanill2004D . The source is also de- 
tected in the first 9 months data of Swift Burst Aler t 
Telescope (BAT) survey of AGNs (|Tueller et al.ll2008h . 
Though being the brightest BLRG in the hard X- 
ray sky above 10 keV, 4C 50.55, located at (l,b) = 
(93.32, 0.3937), had been unrecognized as a bright X-ray 
source until the INTEGRAL and Swift era, due to the 
obscuration by the Galactic plane. This source is also of 
great interest for understanding the accretion flow onto 
supermassive black holes at high fractions of Eddington 
lumino sity, w hich is estimated to be Lboi/^Edd ~ 0.4 (see 
section 133)1 . IMassetti et al.l (|2004D determined the red- 
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shift to be z = 0.020 ± 0.001 based on the detection of a 
broad Ha emission line in the optical spectrum. From the 
observed flux densities at 1 . 4 GHz with the Very Large 
Array (VLA). [Molina et all (|2007t) estimate the inclina- 
tion angle to be 8 ~ 35°, assuming a Lorentz factor of 
7 = 5 and a typical power rati o between the core an d 
total known for radio galaxies ()Giovannini "et~aLl[l98l . 
The inferred viewing geometry of the nucleus is well con- 
sistent with the classification of 4C 50.55 as a BLRG. 

The initial results o n the X-ray spectra are reported 
bv lMolina et al.l (|2007D . who used the XMM-Newton and 
Swift/XRT data combined with a time-averaged INTE- 
GRAL spectrum. They find that the spectra below 10 
keV are apparently hard, and complex, multiple layers 
of absorber are required. Iron-K features are not sig- 
nificantly detected and the strength of reflection com- 
ponents, R = n/2Tr, where il is the solid angle of the 
reflector, is not tightly constrained (R ;$ 0.9). In this pa- 
per, we present the first simultaneous bro ad band X-ray 
data of 4C 50.55 observed with Suzaku (Mi tsuda et al.l 
I2007D . obtained in a different epoch of the XMM-Newton 
observation. Suzaku carries (then available) three CCD 
cameras called the X-ray Imaging Spectrometers, XIS-0, 
XIS-3 (Front-side Illuminated XIS; XIS-FI), and XIS-1 
(Back-side Illuminated XIS; XIS-BI), and a collimated- 
type instrument called the hard X-ray detector (HXD), 
which consists of Si PIN photo-diodes and GSO scin- 
tillation counters. The XISs and HXD-PIN covers the 
0.2-12 keV and 10-70 keV, respectively. Since AGNs 
are time variable, the simultaneous coverage is critical 
for studying the continuum shape over the broad band, 
in particular to accurately constrain the Compton reflec- 
tion component. Moreover, the Suzaku exposure is quite 
deep (a net exposure of ~ 100 ksec), and thus provides 
us with the best quality dataset so far obtained from this 
source. 

The organization of our paper is as follows. First, we 
describe our Suzaku observation and data reduction in 
section 2. Next, we present the analysis and results in 
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section 3. The results of the Swift/BAT data are also 
presented, whose spectrum averaged over 22 months is 
utilized to constrain the highest energy band up to 200 
keV. We plot the spectral energy distribution of 4C 50.55 
from the radio to Gamma-ray bands to discuss the pos- 
sible contribution of the jet component. Finally, we dis- 
cuss our results in comparison with the previous studies 
of this target and other radio galaxies in section 4. In all 
spectral analysis, we apply the Galactic absorption fixed 
at N^ al = 1.0 x 10 22 cm" 2 (|Kalberla et all 127)051 ) . The 
cosmological para meters (Hp, fl m , Q\) = (71 km s -1 
Mpc" 1 , 0.27, 0.73: iKomatsu et al J 120091 ) are adopted in 
calculating the luminosities. The errors attached to spec- 
tral parameters correspond to those at 90% confidence 
limits. 

2. OBSERVATIONS AND DATA REDUCTION 

4C 50.55 was observed with Suzaku from 2007 April 
16 14:05 (UT) to April 18 11:04 (UT) (observation ID 
702027010), focused on the nominal center position of 
the Hard X-ray Detector (HXD). The net exposures after 
data screening, described below, are 85.0 ksec (XIS-0, 1), 
80.0 ksec (XIS-3), 54.7 ksec (PIN). For the HXD, we use 
only data collected with the PIN diodes, because we find 
that the signal from the source in the GSO data is not 
significant over the systema tic error of » 2% in th e non 
X-ray background (NXB) (jFukazawa et all 120081 ). We 
use FTOOLS (heasoft version 6.6.2) to extract data, and 
XSPEC version 11.3.2ag for spectral fitting. 

2.1. XIS Data Reduction 

To apply the latest calibration, we reprocess the unfil- 
tered event files of the XIS data according to the pro- 
cedures written in The Suzaku Data Reduction Guide 
(ABC Guide). We select the data where the time since 
the South Atlantic Anomaly (T_SAA) passage is longer 
than 436 sec, the elevation angle (ELV) is larger than 5°, 
and the dye-elevation angle (DYE-ELV) is larger than 
20°. The XIS events are extracted from a circular re- 
gion centered on the source peak, and the background 
is taken from a source-free region with the same dis- 
tance from the optical axis as the target. We generate 
RMF files with xisrmfqen , and ARF files with xissimar- 
fgen (|Ishisaki et al.l l2007h . We examine the spectra of 
the 55 Fe calibration source (producing an Mn Ka line 
at 5.895 keV) located on the corners of the XIS chips 
to check the accuracy of the energy response. By fit- 
ting them with Gaussians, we are able to reproduce the 
right central energy within the statistical errors and line 
widths of 3.9 (< 19.2) eV (XIS-0), 0.0 (< 15.6) eV (XIS- 
1), and 0.0 (< 14.0) eV (XIS-3). This verifies that both 
energy scale and energy resolution in the responses are 
well reliable. 

2.2. HXD-PIN Data Reduction 

We also reprocess the unfiltered event files of the HXD 
data as well. This include the time assignment (with 
hxdtime), gain correction (hxdpi), and grade classifica- 
tion (hxdgrade). We select data where the time after 
the SAA passage is longer than 500 sec, a cutoff rigid- 
ity (COR) is larger than 6 GV, and the elevation from 
the earth is larger than 5°. We utilize the "tuned" NXB 
event files provided by the HXD team to produce the 
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Fig. 1. — Long term light curve of 4C 50.55 in the 14-195 keV 
band obtained with Swift/BAT with 16 days bin. The fluxes are 
converted from the count rate by assuming a power law photon 
index of V = 1.7. The dotted lines indicate the period of the 
Suzaku observation (from 2007 April 16 to April 18). 

background spectra, to which the cosmi c X-ray back- 
groun d (CXB) based on the formula by iGruber et al.1 
( 1999 1) is added. We use the HXD/PIN response file 
ae_hxd_pinhxnome3_20080129.rsp. 

3. ANALYSIS AND RESULTS 

3.1. Swift/BAT Data 

Figure [1] shows the long term light curve of 4C 50.55 
in the 14-195 keV band obtained with Swift/BAT from 
2004 December 15 to 2009 October 10. Each data point 
corresponds to the averaged flux for 16 days. Time 
variability is clearly noticed. The epoch of our Suzaku 
observation is indicated by the dashed line in the fig- 
ure. The time-averaged BAT spectrum over the first 22 
months covering the 14-195 keV band is plotted in Fig- 
ure [H We find that it can be fit with a single power 
law of T = 1.68 ± 0.25. The time averaged 14-195 keV 
flux is 1.7 x 10~ 10 ergs cm~ 2 s _1 (based on the best-fit 
power law model), which corresponds to a luminosity of 
1.5 x 10 44 ergs s _1 in the rest-frame 14.3-199 keV band. 
In the subsequent subsections, we will apply more real- 
istic spectral models to this spectrum. 

3.2. Suzaku Light Curve 

We make the Suzaku light curves of 4C 50.55 with 
a 5760 sec bin, the orbital period of the satellite, to 
remove any possible modulation related to the orbital 
condition. Figure [3] shows the X-ray light curves ob- 
tained with the XIS-FIs (2-10 keV, upper) and with the 
HXD/PIN (15-40 keV, middle), and their hardness ratio 
(PIN/XIS, lower). The zero point of time corresponds 
to the start time of the Suzaku observation. The XIS 
light curve in the 2-10 keV band indicates a rapid flux 
increase by ~ 20% in the last 20 ksec exposure. The time 
scale of this variability is ~ 10 4 sec, indicating that the 
emission region is within ~ 30r g (r g = is the gravita- 
tional radius) for a black hole mass of M = 10 78 M Q (see 
section |3~51) . By contrast, the light curve in the 15-40 
keV band does not show evidence for significant variabil- 
ity above the statistical errors in the same epoch. Thus, 
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Fig. 2. — The time-averaged Swift/BAT spectrum of 4C 50.55 
over 22 months. The best-fit model is a pexriv model (see Table[T]l. 
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Fig. 3. — Light curves of 4C 50.55 obtained by the Suzaku 
observation, with a 5760 sec bin. The start time corresponds to 
2007 April 16. The dotted line shows the border between epochs 
1 and 2. (upper) the light curve of XIS-FIs in the 2—10 keV 
band, (middle) that of PIN in the 15 - 40 keV band, (lower) the 
hardness ratio between them (15-40 keV / 2-10 keV). 

there is a hint for decrease in the hardness ratio between 
the PIN and XIS, although its significance is marginal. 
For the following spectral analysis, we separate the ob- 
servation period into two, epoch 1 (0-1.4xl0 5 sec) and 
epoch 2 (1.4xl0 5 -1.7xl0 5 sec). 

3.3. Spectral Analysis with Phenomenological Models 
3.3.1. Individual Fit to the Suzaku and Swift Data 

We analyze the Suzaku spectra of epoch 1 and epoch 2 
separately, by performing simultaneous fit to those of the 
XIS-FIs in the 1-9 keV band, the XIS-BI in the 1-8 keV 
band, and the HXD/PIN in the 12-60 keV band. Given 
the fact that the spectrum of 4C 50.55 is subject to a 
heavy Galactic absorption, we do not utilize the XIS-BI 
data below 1 keV to avoid any possible uncertainties in 
the response. We find that the effects of pile-up are sig- 
nificant above 9.0 keV for the XIS data. Further, we 
discard the XIS data of the 1.7-1.9 keV band because of 
calibration uncertainties associated with the instrumen- 



tal Si-K edge. In the simultaneous fit, the flux normal- 
izations for the XIS-FIs and XIS-BI are set free, while 
that of t he PIN is fixed at 1.18 relative to that of the 
XIS-FIs (IMaeda et al.l[2008h . 
We firstly apply a single power law model (modified by 

the Galactic absorption fixed at N^" 1 — 1.0 x 10 22 cm -2 ) 
to the XIS and PIN data in epoch 1, covering the 1- 
60 keV band. The fit is found to be far from accept- 
able (x 2 / dof = 13035/810) and results in a very flat 
slope, r w 1.2. This is because the spectral shape be- 
low 10 keV is much harder than that above 10 keV. 
Next we apply a cut-off power law model, in the form 
of E~ r x exp(— E/E cu t), over which a narrow Gaus- 
sian is added to represent an iron-K emission line at 
6.4 keV. We obtain T « 0.83 and E cut « 14 keV with 
X 2 /dof = 4616/808. There remain strong absorption 
features around 1 keV suggesting intrinsic absorption, 
however. When we add another absorption at the source 
rcdshift (z = 0.02), the fit becomes much better, yield- 
ing T = 1.39 ± 0.02, E cut = 45 ± 2 keV, and N H = 
(0.67±0.01) x 10 22 cm" 2 with y 2 /dof = 997/80 7. Finally, 
similar to the analysis done by I Molina et al.l (|2007t) , we 
consider double absorber model with two different col- 
umn densities, and N^, whose covering fraction is / 
and (1 — /), respectively. The fit is significantly improved 
( X 2 /dof = 940/805) with V = 1.61+°/° 3 5 , E cut = UOtfo 

keV, 7V^« 8.3 x 10 22 cm~ 2 , 7V 2 « 0.73 x 10 22 cm" 2 , and 
/ = 0.19. The photon index becomes a reasonable value 
for AGNs. Note that our "double absorber" model has 
three free parameters for the absorber, while the "dou- 
ble partial covering" mod el ("pcf*pcf*" in th e XSPEC 
terminology) adopted by IMolina et al.l (|2007|) has four. 
Since the latter model does not give a significant im- 
provement for our data (Ax 2 ~ 1), we adopt the former 
as a base continuum model for the following analysis. 

The high quality Suzaku spectra are quite useful to 
constrain the reflection component, which is indicated 
by the presence of the iron-K emission line. Thus, we 
include it by uti lizing our modified version o f "pexriv" 
reflection code (Ma gdziarz & Zdziarski 1995) that as- 
sumes a cutoff power law continuum and contains a self- 
consistent fluorescence iron-K line c alculated accordin g 
to the same algorithm as described in lZvcki et al.l (1999). 
The additional free parameter is the relative reflection 
strength, R(= 0/2tt), while we fix the ionization param- 
eter, temper ature, and inclinatio n angle at 0, 10 5 K, and 
35 degrees (jMolina et al.ll2007f). respec t ively. The so- 
lar Fe abundance by I Anders k, Grevessd (|1989l ) (Fe/H = 
4.68 x 10~ 5 ) is assumed. Considering that the reflec- 
tion most likely occurs in the accretion disk, we smear 
both reflected conti nuum and iron-K e mission line by the 
"diskline" kernel (|Fabian et al.l [19891) . The innermost 
radius r; n is set as a free parameter by assuming an emis- 
sivity law of r~ 3 for a fixed outer radius r out = 10 5 r g . 
When constraining r- m , we utilize only the XIS spectra 
around the iron-K band (3-9 keV for the XIS-FIs and 
3-8 keV for the XIS-BI) and fix all the other parameters 
except for the normalization. Thus, the spectral fit is 
performed by iteration; after rj n is determined from the 
XlS-only fit, it is then fixed when finally determining the 
continuum parameters in the XIS+PIN fit. 

The results of the spectral fit to the individual Suzaku 
spectra in epochs 1 and 2 are summarized in Table [TJ 
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TABLE 1 

Best-Fit Parameters of Cut-off Power Law Model to 
the Individual Suzaku and Swift Data 



TABLE 2 

Best-Fit Parameters of Cut-off Power Law Model to the 
Combined Suzaku and Swift Data 



Parameters 


Epoch 1 


Epoch 2 


Swift /BAT 


Parameters 


Epoch 1 


Epoch 2 (1) 


Epoch 2 (2) 


(a) f 

72-10 kcV 


*7 /t w in — 11 
I.4X 1U 


Q c; w in — 11 
O.O X 1U 




(a) r 

1 '-^2-10 kcV 


o. 1 


("1 Q 

y.o 


y.u 


(a) f 

J 10 — 60 kcV 

^/l4-195 kcV 


1.6 x 10 _lt) 


1.6 x 10 — 10 


1.7 x 10~ 10 


r 


1.0 § 
1.65 ±0.04 


1.0 s 
1.69 ±0.02 


1 n§ 
±.u° 

1.65 § 


(b^gal 

r 

Bcut [keV] 
R{= Q/27r) 
(c) / 

EW [eV] 


1.0 § 
1.61 ±0.05 

80±f« 
0.18 ±0.04 
19+ 003 
76 +i.5 
'■ D -1.4 


1.0 § 
1.64 ±0.03 

* 

0.19§ 
7.6§ 


i fio+0.23 
i - b8 -0.25 

127tlf 
0.18 § 


Scut [keV] 
i?(= Q/2tt) 
_EW [eV] (Fe Ka line) 

(c) / 

(d)7v2 


105 +2 J 
0.17 ±0.04 
21 

0.21 ±0.03 

7 9+1-4 
' - y -1.3 

0.75 ±0.03 


-i Qy + 25 

iz '-23 
* 

18 
0.21§ 
7.9 § 
0.75 § 


105§ 
* 

18 

0.40±°;?° 

2 9+ 16 
z - a -i.o 

49+ ' 17 


0.73 ± 0.03 
22 


0.73 s 
19 




*7 dof 


926.9/812 


201.1/182 


192.7/181 


r-m[rg] 


720§ (> 340) 1 " 


720 § 


720 § 


Note. — Errors are 90% confidence level for a sinj 


;le parameter. 



916.0/805 



190.8/175 8.594/5 



Note. — Errors are 90% confidence level for a single 
parameter. 

( a ) Observed fluxes in the 2-10 keV, 10-60 keV, and 14-195 
keV bands, in units of ergs cm -2 s — 1 . 

( b ) Galactic absorption column density in units of 10 22 cm -2 . 

( c ) Covering fraction 

( d ) Intrinsic absorption column density at the source redshift 
in units of 10 22 cm -2 . 

s Parameters fixed at these values. 

* We assume the same reflection component as that deter- 
mined in epoch 1. 

t We constrain the inner radius only from the 3 — 9 keV XIS 
spectra in epoch 1, which is then fixed at the best-fit when 
determining the continuum parameters. 

Figure |4] shows the XIS+PIN spectra folded by the en- 
ergy responses, over which the best-fit models are plot- 
ted, with residuals in the lower panel. The expanded 
figure of the XIS spectra between 3-9 keV in epoch 1 is 
plotted in Figure [SJ To emphasize the iron-K line fea- 
ture, the residuals when the line is excluded from the 
model are shown in the lower panel of this figure. From 
epoch 1, we obtain V = 1.61 ± 0.05, E cut = 80±f| keV, 
and i? = 0.18±0.04. The innermost radius is constrained 
to be r in = 720r g (> 340r g ) from the XIS data. For the 
analysis of the epoch 2 spectra, we assume the same pa- 
rameters of the reflection component (including its abso- 
lute flux) as those found from the epoch 1 data, since it is 
very unlikely that it varied on such a short time scale of 
< 10 4 sec. The parameters of the absorption are fixed to 
the epoch 1 values as well. Thus, only free parameters are 
r, Scut; an d the normalization. Finally, we also perform 
spectral fit to the time-averaged Swift/BAT spectrum by 
adopting the same model. The reflection strength is fixed 
at R = 0.18 referring to the epoch 1 result. The best-fit 
model is over-plotted in Figure [2 whose parameters are 
summarized in Table [1] 



3.3.2. Simultaneous Fit to the Suzaku and Swift Data 

From the above analysis, we find no significant differ- 
ences in the spectral parameters (except for the normal- 
ization) within the statistical errors between the Suzaku 
epoch 1, epoch 2, and Swift/BAT data, although there 
is a hint that the spectrum became slightly softer in 
epoch 2. Thus, to best constrain the continuum parame- 
ters, in particular the cutoff energy, we study the Suzaku 



( a ) Intrinsic luminosity in the 2-10 keV band corrected for both 
Galactic and intrinsic absorptions in units of 10 43 ergs s" 1 . 
Galactic absorption column density in units of 10 22 cm" 2 . 

( c ) Covering fraction. 

( d ) Intrinsic absorption column density in units of 10 22 cm" 2 . 
§ Parameters fixed at these values. 

* We assume the same reflection component as that determined in 
epoch 1. 

spectra (either of the two epochs) in the 1-60 keV band 
and Swift spectrum in the 14-195 keV band simultane- 
ously, in all following analysis. The flux normalization 
between the Suzaku (XIS-FIs) and BAT spectra are set 
free, to take into account the time variability. In analyz- 
ing the spectra of epoch 2, we always fix the reflection 
component to that determined from the epoch 1 data. 

Table [2] summarizes the results using the same phe- 
nomenological model (cutoff power law) as adopted in 
section l3~5l For epoch 2, we consider two extreme cases 
as the cause of the spectral variability from epoch 1 that 
(1) only the continuum changed without change of the 
absorber and that (2) only the absorber changed with the 
same continuum except for its normalization. We obtain 
similarly good fits for the two cases, and thus both pos- 
sibilities are plausible from the spectral analysis. In re- 
ality, however, it may be difficult to explain such a short 
time (< 10 4 sec) variability by the absorber alone. If the 
absorber makes Kepler motion at ^1000 r g , a typical lo- 
cation of the broad line region in AGNs, it moves only ~ 
r g in 10 4 sec, by assuming the black hole mass of 10 7 8 
Mq. Thus, unless the emitting region is extremely small 
(like < several r g ), it is unlikely that crossing blobs in the 
line of sight can cause the large variability as observed. 

3.4. Spectral Analysis with Comptonization Model 

In this subsection, we analyze the spectra of 4C 50.55 
with a physically motivated model instead of the phe- 
nomenological "cutoff power law" model, which is a 
mathematical approximation of the X-ray spectra of 
AGNs. Such analysis of AGN spectra has been very 
limited so far, since it requires high quality broad band 
data. As we will discuss in section T3.5[ the contribution 
from the jet components is very small in the X-ray band. 
Hence, we consider that the origin of the continuum emis- 
sion is predominantly thermal Comptonization of soft 
(ultra-violet) photons off hot electrons in the corona lo- 
cated above the accretion disk. Accordingly, we adopt a 
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Fig. 4. — (a) (left) The folded spectra of 4C 50.55 in epoch 1 obtained with the XIS-FIs in the 1-9 keV band (black), the XIS-BI in the 
1-8 keV band (red, open circle), and the PIN in the 12-60 keV band (blue). The solid curve represents the best-fit model (see Table [TJ- 
The dotted curve represents the reflection component, (b) (right) those in epoch 2. 
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Fig. 5. — (upper) The spectra of 4C 50.55 obtained with the 
XIS-FIs (black) and XIS-BI (red, open circle) in epoch 1. The 
best-fit model (see Table [TJ is plotted by the solid curve. The 
dashed curve indicates the model from which the iron-K emission 
line is excluded, (lower) The residuals of the fit in units of \ for 
the model without the iron-K line. 



therm al Comptonization model, thComp (Zvc kTet al.l 
1999), for the primary continuum. It has two free pa- 
rameters, the slope T and electron temperature kT c . The 
electron scattering optical depth r e is related to T e and T 
by the following formula (jSunvaev fc Titarchukl ll980): 



'2.25 



(Te/511 keV)[(r + 0.5) 2 - 2.25] 



1.5(1) 



For seed photons, we assume a multicolor disk compo- 
nent with the innermost temperature of 0.01 keV. This 
choice is not important to constrain the Comptoniza- 
tion parameters. As described in the previous subsec- 
tion, the reflection component and fluorescence iron-K 
line are self-consistently included, which are blurred by 
the "diskline" profile with the same parameters as ob- 
tained above (r; n = 720r g ). 

We perform simultaneous fit to the Suzaku (epoch 1 or 
2) and Swift/BAT spectra using this model. The best-fit 



TABLE 3 

Best-Fit Parameters of Comptonization Model to the 
Combined Suzaku and Swift Data 



Parameters 


Epoch 1 & BAT 


Epoch 2 & BAT 


(a) T 

k '^2-10 keV 


8.4 


9.6 


(b )Ar gal 


1.0 § 


1.0 § 


r 


1.78 ± 0.02 


1.81 ± 0.01 


kT a [keV] 


31±J° 




<°>( r c 




2 1+ ' 5 1 
z - x -0.8 > 


R(= U/2tt) 


0.16 ±0.04 


* 


EW [eV] (Fe Ka line) 


19 


17 


(d)/ 


0.27 ±0.02 


0.27§ 






8.6§ 




0.80 ±0.03 


0.80§ 


X 2 / dof 


942.5/812 


200.8/182 



Note. — Errors are 90% confidence level for a single 
parameter. 

( a ) Intrinsic luminosity in the 2—10 keV band corrected 
for both Galactic and intrinsic absorptions in units of 
10 43 ergs s — 1 . 

' b ) Galactic absorption column density in units of 10 22 cm" 2 . 
' c ) Electron-scattering optical depth calculated from the 
equation (1). 
( d ) Covering fraction. 

( e ' Intrinsic absorption column density in units of 10 22 cm" 2 . 
§ Parameters fixed at these values. 

* We assume the same reflection component as that deter- 
mined in epoch 1. 

parameters are summarized in Table [3] We obtain the 
electron temperature of «30 keV with an optical depth 
of « 3. Again, the reflection component in epoch 2 are 
fixed at the same one determined from epoch 1. We find 
that the continuum parameters are consistent each other 
between epochs 1 and 2 within the statistical errors at 
the 90% confidence level. Figure [6] shows the results in 
the two epochs, in the form of unfolded spectra (i.e.,, 
those corrected for the effective area of the instruments) 
in units of EI(E), where E is photon energy and 1(E) 
the energy flux. 

3.5. Multi- Wavelengths Spectral Energy Distribution 
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Fig. 6. — (a) (left) The unfolded spectra of 4C 50.55 determined from the Suzaku data in epoch 1 and from the Swift/HAT data, based on 
a thermal Comptonization model. The crosses (black) represent the data points, dashed curve (blue) the transmitted component, dashed 
curve (magenta) the reflection component, and solid curve (red) the total. The best-fit parameters are given in Table |3] (b) (right) those 
in epoch 2. 

sion in 4C 50.55 is simila r to that of jet dominated- 
sources (jKubo et alJ 119981 ) with correction for an esti- 
mated beaming factor (5 ~ 1 for 4C 50.55 and 8 ~ 10 
for blazars), the predicted X-ray emission is a factor of 
~ 10 1-4 less than the observed X-ray flux. We thus con- 
clude that X-ray emission due to a jet via synchrotron 
emission or inverse Compton should be very small in 
the total X-ray emission of this source. Based on the 
Suzaku results, the ratio of the core luminosity at 5 GHz 
and that in the 2-10 keV band is estimated to be log 
i?x = —3.6. Thus, while 4C 50.55 should be cl assified as 
a radio loud object (jTerashima fc Wil son 2003J), its radio 
to X-ray power ratio is much lower compared with typi- 
cal blazars and mo re powerful radio gala xies, like 3C 120 
(log R x = -2.1) dKataoka et all 120071) . This confirms 
the conclusion bv lMolina et al.l ( 20071 ). who did not find 
any evidence for an additional power law component rep- 
resenting the jet emission in the XMM-Newton spectrum. 

We note that the optical fluxes of 4C 50.55 are highly 
uncertain at present due to the extinction (both Galac- 
tic and intrinsic ones), which could largely affect the 
empirical estimate of the black hole mass. In fact, 
the extinct ion-corrected AGN continuum flux at 5100 A 
adopted bv lWinter et al.1 (|2010D falls far below an extrap- 
olation from the X-ray spectra toward lower energies. 
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Fig. 7.— The SED of 4C 50.55. The solid (black) and dotted 
(red) curves represent the Suzaku spectra in epochs 1 and 2, re- 
spectively, while the thick (purple) one represents their difference 
spectrum between the two epochs. The spectrum of Sun/i/BAT is 
shown by the blue crosses. The pink arrow re presents the uppe r 
limit on the 7-ray flux from th e Fermi data I Abdo et al. 2009). 
The radio data are taken after Molina et al. I ( 120071) ; the orange 
filled-triangle and reddish v iolet filled-circles rep resent the core flux 
obtained from the GMRT QPandev et al.112009 1 and the VLA, re- 
spectively. The blue filled-square represent the total flux including 
both core and lobes. 



We summarize in Figure [7] the available multi- 
wavelengths fluxes of 4C 50.55 to discuss its spectral en- 
ergy distribution (SED) from the radio to 7-ray bands. 
In the radio band, we utilize those obtained by the VLA 
and the Giant Metrewave Rad io Telescope (GMRT) com- 
piled by IMolina et al.l (|2007l) . In the GeV 7-ray band, 
the upper limit d erived from the one-year Fermi data 
(|Abdo et all I2009D is indicated by the arrow. We plot 
the results of cutoff power law fit to the Suzaku spec- 
tra in epochs 1 (black) and 2 (red) as well as the time- 
averaged Swift/BAT data, correcting for both Galactic 
and intrinsic absorptions. To examine the origin of the 
time variability, we analyze the difference spectrum of 
XIS-FIs between epochs 2 and 1. Fitting with a single 
power law with the same double absorption as given in 
Tabled! we find T « 1.8. This result is also plotted in 
the figure (green curve). 

Assuming that the global intrinsic SED of the jet emis- 



As a more simple approach, assuming typ ical SEDs of 
AGNs ijElvis et alJliOM IGnme et alJl201f)l) . we roughly 
estimate the 5100 A flux to be (0.7 - 8) x 10" 10 erg 
cm~ 2 s~ 4 from the BAT flux averaged for 22 months. 
Then, from the luminosity at 5100 A, (0.6 - 7) x 10 44 
erg s _1 , and the width of the line, 2320 km s _1 
FWHM (|Winter et al.l 120101) . we derive the black hole 
mass of 4C 50. 55 to be ~ iq7-8±Q-3 m usin g t ^ e f or _ 
mula given by iVestergaard fc Peterson! (|2006h . This is 
ab out 10 times higher than the estimate (10 6 - 58±0 07 M ) 
bv lWinter et al.ll2010L Nevertheless, the suggestion that 
the black hole in 4C 50.55 accretes with a high Edding- 
ton ratio still holds. With a typical bolomctric correction 
factor (iboi/^2-10 keV =30), we estimate the fraction of 
Eddington luminosity of 4C 50.55 to be Lboi/^Edd ~ 0.4. 

4. SUMMARY AND DISCUSSION 



Suzaku Observation of 4C 50.55 
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We have obtained the first simultaneous, broad band 
X-ray spectra of 4C 50.55 over the 1-60 keV band with 
Suzaku. These provide one of the best quality high en- 
ergy data so far obtained from BLRGs. From the SED, 
we conclude that there is little contribution from jets in 
the X-ray emission, which is thus dominated by Comp- 
tonization by hot corona. We find that the overall contin- 
uum is represented by a cutoff power law, or a thermal 
Comptonization model, with complex intrinsic absorp- 
tions. We show that at least two absorber with different 
covering fractions and column densities are necessary to 
model the spectrum. The photon index and cutoff en- 
ergy of 4C 50.55 are wit hin the distribution of previous 
observations of BLRGs (|Grandi et al.l 1200611 , although 
both are somewhat smaller (r ss 1.6 and E cut ~ 100 
keV) compared with typical values foun d from Seyfer t 
galaxies (r wl.9 and E cut ~200 keV, iDadinal 120081 ). 
These find ings on 4C 50 . 55 are fully consistent with the 
results by Molina et all (|2007l ) from the XMM-Newton 
and INTEGRAL data observed on earlier epochs. We 
estimate that the corona responsible for Comptoniza- 
tion is optically thick for scattering r e w 3 and is rel- 
atively cool T c pa 30 keV. The optical depth is larger 
and the temperature is lower than those obtained from 
Seyfe rt galaxies with ther mal Comptonization models 
(e.g., iLubinski et al.|[20Tol) . as expected from the com- 
parison in r and E cut - 

From our observations, time variability of the X-ray 
flux on both long ( £ 10 6 sec) and short (~ 10 4 sec) 
time scales is detected. The averaged 2-10 keV flux was 
the highest in our Suzaku observa tion among those in the 
previous observations reported in lMolina et al.l (|2007l ) by 
a factor of 1.3-2.5. The hard X-ray flux in the 17-100 
keV band (1.8 x 10~ 10 erg cm~ 2 s" 1 ) was also higher by 
1.4 than the averaged flux over 22-months obtained with 
Swift/BAT (1.3 x lO" 10 erg cm" 2 s" 1 ). The large long- 
term variability in the hard X-ray above 10 keV suggests 
that it is mainly produced by the intrinsic emission, not 
purely by the chan ge of the absorber as discussed by 
iMolina et al.1 (l2007fl in the frame w ork of a "patch torus" 
model (|Elitzur fc Shlosmanll2006[ ). While there remains 
a possibility that Compton thick absorber (Ah > 10 24 
cm -2 ) may completely cover a part of the emission region 
to cause the flux variability, it would produce signals of 
heavy obscuration such as a deep iron-K edge feature in 
the X-ray spectra, which are not seen in the data. Thus, 
at least variability of the continuum flux is required to 
explain these results, while that of the absorber can also 
contribute to t he varia bility , in particular below 10 keV 
(|Risaliti et al.l see e.g., [20051 for NGC 1365). 

We significantly detect an iron-K emission line and 
obtain a tight constraint on the reflectio n component, 
even t hough it is quite weak as reported bv IMolina et al.l 
(2007). The reflection strength, R ~ 0.2 is much smaller 
compar ed with Seyfe rt 1 galaxies, which typically have 
R ~ 1 (Dadina 2008). Even correcting for a variability 
effect that the direct continuum flux was higher in our 
observations by a factor of «1.4 than the 22-months aver- 
age, the reflection is still small, R ~ 0.3. The weak reflec- 
tion is also consistent with the narrow iron-K emission 
line, which indicates that the reflection is mainly pro- 
duced by relatively outer parts of the disk (hence with a 
small solid angle) , unlike the results from typical Seyfert 



1 galaxies (Dadina 200 8]). Our 4C 50.5 5 resul t confirms 
the trend reported bv iSambruna et ail ([2002) for radio 
loud AGNs. 

The analysis of the iron-K line profile yields an (appar- 
ently) large innermost radius, r m ~ 700r g , by assuming 
an emissivity law of r~ 3 . Our diskline result suggests it 
unlikely that a "standard disk" extends down to close 
to the innermost stable circular orbit (ISCO) around the 
black hole, < 6 r g . Instead, it is possible that the disk is 
there but its inner part is covered by an optically thick 
corona, as estimated by our Comptonization model fit, 
which would smear out relativistic iron-K broad lines. 
Note that the obtained r- ln value does not directly mean 
that the disk is truncated at that radius, because (1) the 
estimated r; n critically depends on the emissivity profile 
and (2) there may be another line component from dis- 
tant parts, such as the torus. If the scale height of the 
X-ray irradiating corona is sufficiently small, then one 
would expect a flatter slope for the emissivity law, even 
close to r~ 2 . In this extreme case, we were not able to 
obtain good constraints on r- m from our data. To exam- 
ine the second possibility, we apply a two-components 
line model to the XIS 3-9 keV spectra, consisting of a 
narrow Gaussian at 6.4 keV and a broad diskline, which 
represents that from the torus and disk, respectively. We 
obtain a worse fit than the single diskline fit by A% 2 = 5 
even with a larger degrees of freedom, suggesting that the 
two components model is not a good description of the 
data. Nevertheless, when r- 1B of the disk l ine component 
is fix ed at 10 r g as found from 3C 120 (Kataoka et al. 
l200l . both lines are found to be significant with equiv- 
alent widths of 19±6 eV and 24±15 eV, respectively. 
Thus, we do not completely exclude the possibility for 
the presence of a moderately broadened iron-K line in 
the observed spectra of 4C 50.55. 

The inferred geometry of the accretion disk in 4C 50.55 
(i.e., truncated and/or inner parts covered by corona) 
may be common features of AGNs with powerful jets. 
Recent Suzaku studies indicate that radio galaxies also 
have relatively narr ow iron-K emission lin es e.g., r; n > 
20 r g for 3C 390.3 (ISambruna et all 120091) and r in > 44 
r g for 4C +74.26 (jLarsson et all 120081 ) from the single 
diskline fit, and r m = (9 ± l)r g for 3C 120 from the mul- 
tiple components fit (|Kataoka et al.ll2007f ). This result 
is in accordance with an expectation from theories that 
jets are more easily produced by radiatively inefficient 
accretion flow than by a standard disk. 

Another key parameter to understand the accretion 
flow is the Eddington ratio, which is estimated to be 
^boi/^Edd ~ 0.4 for 4C 50.55 (section |3~5"| . Similarly, 
we also estimate that of 3 C 120 to be Lh^ /L^ dd ~ 0.5, 
using the 2-10 keV flux (jKataoka et alJl2007Q and the 
black hole mass of 10 7 7 M (|Peterson et al.ll2004D . Thus, 
these two sources may belong to a very similar class of 
AGNs, except for the radio loudness to the X-ray flux (log 
Rx = -2.1 for 3C 120 and log R x = -3.6 for 4C 50.55), 
which could be partially explained by the small inclina- 
tion angle of 3C 120 (i < 14°; see lKataoka et all 120071) 
compared with 4C 50.55 (i ~ 35°). The physical reason 
for the difference in their X-ray spectra that the reflec- 
tion component and iron-K lines are stronger in 3C 120 
(R ~ 0.7) is not clear at present. 

4C 50.55 and 3C 120 are rare objects having 
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distinctively high fractions of Eddington luminosity 
compared with other typical BLRG s, for instance, 
L hn] /L FlM = 0.01-0.07 for 3C 390.3 (|Sambruna et al.1 
2009; iLewis fc Eracleousl [2006h. ~ 0.04 for 4C +74.26 
([Larsson et all I2008D. an d 0.001-0.002 for Arp 102B 
([Lewis fc Eracle ous 2006J). By analogy to the Galactic 
black holes, these low Eddington ratio sources likely cor- 
respond to the low/hard state, where the accretion disk is 
accompanied by steady jets, while normal Seyfert galax- 
ies may do to the high/soft state, where the disk extends 
close to the ISCO with quenched jet activity. The ac- 
cretion flows in 4C 50.55 and 3C 120 could be explained 
as a high luminosity end of the low/hard state. Alter- 
natively, they may be another state achieved with even 
higher mass accretion rates than in the high/soft state, 
where the disk structure is also similar to that found 
in the low/hard state (i.e., truncated disk). For this 
possibility, it is interesting to note the similarly to the 
high-Eddington ratio Galactic black hole GRS 1915+105, 
which exhibits a similarly narrow iron-K emission line 
over a Comptonization dominated continuum, imply- 
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ing that the inne r disk is fully covered by a corona 
(jUeda et all 12010ft ; in GRS 1915+105, a compact jet is 
also detected in a steady stat e with a hard spectrum , 
so-called in Class \ ( see e -g-i iFender fc Bellonil [2004) . 
In summary, the unified picture of accretion flows over 
a wide range of black hole mass is far from established. 
Further systematic studies of the accretion disk structure 
of radio loud AGNs at various accretion rates based on 
detailed X-ray spectroscopy and multi- wavelengths data 
are very important to reveal these fundamental problems. 
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